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Abstract: We present a method for solving BPS equations obtained in the collective- field 
approach to matrix models. The method enables us to find BPS solutions and quantum 
excitations around these solutions in the one-matrix model, and in general for the Calogero 
model. These semiclassical solutions correspond to giant gravitons described by matrix 
models obtained in the framework of AdS/CFT correspondence. The two-field model, 
associated with two types of giant gravitons, is investigated. In this duality-based matrix 
model we find the finite form of the n-soliton solution. The singular limit of this solution 
is examined and a realization of open-closed string duality is proposed. 
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1. Introduction 

In the framework of AdS/CFT correspondence, the dynamics of giant gravitons was studied 
using matrix model with the harmonic-oscillator potential [H, [21 [| . The interpretation of 
the matrix eigenvalues as fermions allows a description of gravitational excitations in the 
holographic dual of N = 4 SYM in terms of droplets in the phase space occupied by 
fermions. In particular, it was shown that giant gravitons expanding along AdS§ and 
S 5 could be interpreted as a single excitation high above the Fermi sea, or as a hole in 
the Fermi sea, respectively [Q]. In Ref.[Q] a correspondence between the general fermionic 
droplet and the classical ansatz for the AdS configuration was established. From Ref.Jl| 
it follows that giant gravitons are described by particular combinations of single-trace 
and multi-trace operators, known as Schur polynomials, pictorially represented by Young 
diagrams. For example, the afore-mentioned giants on AdS§ and S 5 are represented by 
two types of Young tableaux: the one with the one-row diagram and the other with the 
one-column diagram |2|. The matrix model with the harmonic-oscillator potential is 
related to the free matrix model via su(l,l) algebra which contains Hamitonians of both 
models as generators. As a consequence, their eigenstates are related via coherent states 
or by time reparametrization ||. Therefore, a detailed discussion of BPS solutions of the 
matrix model without harmonic potential is relevant for giant graviton physics. 

In this paper we analyse the model without the external potential and present a method 
for obtaining BPS solutions. Our starting point is the background-independent model 
defined by the action 

S=\j dtTrM 2 (t), (1.1) 
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where M is the N x N matrix. In Ref.||, the analysis of the states of the free, U(N) 
gauge-invariant matrix model represented by Young diagrams revealed two single-particle 
branches. The one-row diagram represents a single particle-like excitation above the Fermi 
sea with the dispersion law 

1 _ hi 

(1.2) 



u{k) = \{k + k F Y- k -f, 



and the one-column diagram represents a hole-like excitation in the Fermi sea with the 
dispersion law 



kip 



1 



(k - k F f 



(1.3) 
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The dispersion laws (|1.2| ) and (|1.3j ) made possible the construction of effective collective- 
field Lagrangians for each of these branches ||, which are recognized as Lagrangians of 
0(N) pi] and Sp(N) B gauge-invariant matrix models. These matrix models are connected 
by duality || , similarly to the duality between S 5 and AdS& giant gravitons |]l0[ . As noted 
in Ref . |Ll|] ? the 0(N) and Sp(N) gauge-invariant matrix models have their own string 
interpretation as unoriented (super)strings in two dimensions. 

On a singlet subspace the free matrix model ( |1.1| ) reduces to the quantum mechanics 
of the N eigenvalues Xi of the matrix M. Introduction of the invariant measure over 
the matrix configuration space into the wavefunctions produces a prefactor nl<j( x i ~ x j) X 



which specifies generalized statistics [12], controlled by the parametar A. For matrix models 
the parameter A is related to the number of independent matrix degrees of freedom and 
takes the value A = 1/2, 1, 2, in the case of O(N), U(N) or Sp(N) gauge-invariant matrices, 
respectively ||. Finally, the dynamics of the eigenvalues is determined by the Hamiltonian 
of the Calogero type [O] 



tf C M({*};A) = -i£^ + A(A-l)£ 

i 1 i<j 



(1.4) 



Keeping in mind the relevance of the Calogero-type models for different branches of physics, 
for example in the quantum Hall effect [|l4j], spin models [ji5| 1, hydrodinamical models |l6| . 
systems of particles with generalized statistics |]l2"|| , two-dimensional QCD |l7j and black- 
hole physics |l8|], we use a general A throughout the paper. 

In order to study the large- N, continuum limit of the model defined by action (|1 . 1|) , 
we use the collective-field approach, developed in Ref.[[ll|. Introducing the collective fields 
p and 7r 



p(x,t) = ~J ^e tkt p k (t), p k (t) 
[d x ir(x),p(y)} = -id x 5(x - y), 



-Tre~ ikM(t \ 



(1.5) 



the matrix model ( |1.1| ) was recast in the following Hamiltonian: 



#cm([p];A) 



dxp(x)(d x ir(x)) 2 + 
P 



dx p{x)d x - 



x-y 



dxp(x 
A- 1 



(A-l) 



d x p(x) 
p(x) 



2ir\p H (x) 



y=x 



dx d x 5(x — y) 



y=x ' 



(1.6) 
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where p H is the Hilbert transform of p defined by 

p(y) 



p H {x) = ---fdy- 



x-y 

The collective field p{x) can be viewed as a bosonized 'fermionic' wave function of the afore- 
mentioned discrete case. Another approach to bosonization, suited for a finite number of 



fermions in one space dimension was proposed in Ref.[2C]. 

In attempt to go beyond the case of free fermions, the authors of Ref.|2l| start with 
two matrices and treat one of the matrices in the collective-field theory approach, while 
the other is treated in the coherent-state representation. This leads to the eigenvalue equa- 



tions first found in [22], describing angular degrees of freedom of the single-matrix model. 
Another generalization of the free fermion picture was proposed in Ref.[^], leading to the 
connection between multi-charge BPS operators of N = 4 SYM and the super symmetric 
generalizations of the spin Calogero system. These developments motivated us to look 
also into the Lagrangian, introduced in Refs.|24|, with two collective fields: one field 
associated with real symmetric matrix (with 0{N) invariance) and the other field associ- 
ated with quaternionic matrix (with Sp{N) invariance). This two-field model arises from 
the decomposition of the hermitian matrix into the sum of symmetric and antisymmetric 
matrix and can be thought of as a model of interacting giant gravitons. 

The plan of the paper is the following. In section 2 we present a method for solving 
the BPS equation for the one-matrix model. The construction of the conformal field en- 
ables us to reduce the problem to the Riccati differential equation, which we relate to the 
Benjamin-Ono equation previously found in ||, [l6|. Using the Riccati equation we investi- 
gate quantum excitations around the BPS solutions and related dispersion laws. In section 
3 we analyse the two-field model and present the form of the BPS solution for n solitons. 
The singular limit of some of these solutions is examined and a realization of open/closed 
string duality is proposed. In the last section we summarize the main results and discuss 
some open questions. 

2. One-matrix model 

2.1 Riccati equation and boundary fields 

Here we present a method for constructing BPS solutions of the Hamiltonian (O]). The 



terms in the second line in ( L6) are singular counter terms, which do not contribute in 
the leading order in N. Assuming that the field d x ii{x) vanishes, the leading part of the 
collective-field Hamiltonian (|1.6| ) in the 1/N expansion is given by the effective potential 



y eff = 2 / dxp{x) 



A — 1 d x p{x) H 1 

2 p{x) 



(2.1) 



The effective potential can be rewritten as 

Veff = E 0+ l -J dxp(x) l^ 9 -^ + (f + v ) - X^(x) 2 



(2.2) 
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where the additional term in V e jf is defined as 



Vcot(qx/2 + (p) = lim 



sin(gx + 2ip) 



e->0 coshe — cos{qx + 2ip) 



(2.3) 



In (|2.2| ), Eq represents the terms which are to be subtracted because of the addition of 
Vcot(qx/2 + <p) term into the square brackets: 



En 



g(A - If 
4 

^vrA(A - 1) 



dxd x p(x)Vcot ( ^- + ip 



g 2 (A-i) 2 



dxp(x)Vcot 2 + if 



dxp(x)p H (x)'P cot ( + y> 



qx 



(2.4) 



Assuming the compact support [—L/2, L/2], using the normalization condition J dxp{x) 
N and the identity 



(f H g + fg H ) H = f H g H -fg + fono, 



dx 



(2.5) 



and performing partial integration one obtains 



En 



qN(l - A) 



(1 - A)< ? + 2vrA 



g7rA(A - l)L / 2 H 2 
+ ~ A [P \x) ~ P {x) 



q(X - l) 2 /(/x 
+ 1 p(x)7>cot ^ — + 99 



x=—2ip/q 



L/2 



+ 



-L/2 



(2.6) 



At this point, q and 99 are free parameters to be determined by boundary conditions 
such that the last two terms in ( p.6[ ) should vanish and by the condition that Eq should 
be a non- negative constant. The precise choice of q and ip satisfying these conditions 
determines different solutions of the model and is discussed in the next subsection. Here 
we proceed with the analysis of the effective potential ( |2.2| ). With Eq being a constant, the 
contribution of V e j f to the Hamiltonian is minimized by a solution of the integro-differential 
Bogomol'nyi-type equation 



d x p = qVcot ( + ip ) p + 



Att 
A-l 



2PP H - 



Taking the Hilbert transform of Eq.(2/7), we find the equation for p H , 



_ rr / qx \ H Xtt 

dxP = qPcot [y + v) p -ipo- v— [ 



2 H 2 2 
P ~P ~ Po 



We construct the conformal field containing only the positive frequency part of p 



$ = p H +ip 



1 



7T 



dz- 



p{z) 



z — x — le 



(2.7) 



(2.8) 



(2.9) 



and satisfying the Riccati differential equation 

A7r t9 _ ( qx \ , XnpQ 



A 



-$ 2 + qVcot + $ + J¥± - qp . 



(2.10) 
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The physical interpretation of the Riccati differential equation can be found from the 
relation to the Benjamin-Ono equation [16]. Taking into account the Bogomol'nyi limit, 
we evaluate the fields u + and u~ from Ref.[16|: 

u~ = v^vr^, u + = ^i-Vcot (— + cp) (2.11) 
y/\ 2 V 2 V 

and then plugging ( |2.11| ) into the Benjamin-Ono equation, we obtain Eq. ( |2. 10| ) . 

Obviously, if the field p is a solution of Eq, Q2.7| ), then the field $ necessarily satisfies 
Eq.fllTol). The converse is not true in general, and to obtain the field p from the solution 



of the Riccati equation ( 2.10 ), it is sufficient that the condition 

$H(x) = i$(x)+ Po (2.12) 



holds. In this case, the solution of Eq.( |2.7j ) is simply given as p = —i($> — $*)/2. Equa- 
tion ( |2.10|) can be further transformed into the second-order Schrodinger-like differential 
equation by making the substitution $ = (1 — X)d x v / (Xtiv): 



dlv = .Vcot (f + <p) d x v -^(^l + q j v . (2.13) 
2.2 Semiclassical solutions 

The construction of the Riccati equation enables us to obtain the semiclassical static solu- 
tions. Generally, there are two methods for solving the Riccati equation. One method is 
to construct a general solution from the known particular solution, while the other method 



is based on solving the second-order linear differential equation (2.13). Applying both 
methods, we find some interesting solutions, given in Table [I]. 

In the following analysis we differentiate two possibilities: A < 1 and A > 1. First, 



we discuss the case A < 1. The particular solution of the Riccati equation ( 2.10 ) for non- 
vanishing po is given in the first row of Table ffl. The parameter of this solution satisfies 
the relation 



l + (2.14) 



Taking into account the boundary conditions 

P H {~— ) = P(- — ) = 0, 7>cot(^ + <p) = 0, (2.15) 
q q 4 

we find 

q = 2ttM/L, M e N, (2.16) 

where the number M can be interpreted as the number of solitons. In order to have odd M, 
we take tp = 0, whereas for even M we take ip = tt/2. Taking into account the normalization 
condition, we find 

From the M-soliton solution in the limit L — > oo, keeping po fixed and defining 

b=(l-X)/(Xnp Q ), (2.18) 
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A $ fl (x) p a (x) E 

TT1 SEES l-e'(^+ 2 ^) g(l-A)coth(t/2) l-cos(g:r+2^)~ (1-AK r , ,,,,,,,, 
A< - L A7r(e'-1) l_ e -t e i(?^+2 v ) 2ttA cosh t-cos(< ? x+2 V ) 2£ 5 [AiV iW+^-AJ-NiWZ j 

1-A ix 1-A j 2 (1-A) 3 

Xirb x+ib Xjrb x 2 +b 2 2Ab 2 

x r i ifc(A-l) i+e-y^ fc(A-l) sinhf o 

A ^ 1 2?rA i_ e -i e »fci 2vrA cosh i-cos fcz U 



1-A 1 A-l 

Xtt x+ib A-7T x 2 +b 2 







A ^ 1 ip Po 

Table 1: BPS solutions 



we find the one-soliton solution (M = 1, ip = 0) obtained in Refs.|25], 26]. The energy is 
just the energy of one soliton obtained by taking the corresponding limit. The uniform 
zero-energy solution p(x) = po is obtained in the limit q — > 0, taking ip = n/2. 

Next we discuss the case A > 1. We take q = 0, 92 = vr/2, thus eliminating the term 
■Pcot from Eq.( 2.1C| ), and obtain the general solution of the Riccati equation (2.10) for 



Po = (A 27r j^ , k = 2tvM/L. (2.19) 

It is given in the third row of Table |l], where t is a non-negative free parameter. An 
additional solution could be obtained from the general solution in the limit L — > 00, taking 
t = 2irb/L. In the case t — ► 00, we obtain the constant density solution p(x) = po. Taking 
into account the normalization condition we obtain that the number of solitons M exceeds 
the number of particles ./V giving us the relation 

A = Mj(M - N). (2.20) 

Solitons on the compact support from Table [l] are of the same shape as solitons in 



the Sutherland model [27], thus reflecting the fact that the two models are interrelated 
via the periodicity condition. Using the condensed-matter language, we interpret the M- 
soliton solutions as soliton trains. In the large-M limit, these solutions can be viewed 
as crystal-like structures with periodicity 2ir/q or 2ixjk. The one-soliton solutions are 
regarded as composite particles (quasi-particles). Resembling the holes (lumps), these one- 
soliton solutions enable us to identify the M-soliton solutions found for A < 1 (A > 1) 
as dark (bright) solitons. In Ref. f26[| it was shown that adding the term (1 — A)/(x — z) 
into the effective potential was equivalent to the extraction of the prefactor ni( x i — ^) 1A 
from the wave function of the Hamiltonian ( [01 ). This equivalence enables us to associate 
a quasi-particle located at z with the prefactor of the wave function. Consequently, the 
additional term Vcot(qx/2 + tp) is associated with the prefactor describing M equidistant 
quasi-particles. 

The soliton solutions we have found in the collective-field formulation of the free ma- 



trix model correspond to the particle and hole states in the system (1.4) of nonrelativistic 
fermions ||. Owing to the su(l, 1) dynamical symmetry H, the eigenstates of the Hamilto- 
nian fll.4| ) can be represented as generalized coherent states of the same Hamiltonian with 
the additional harmonic potential interaction between fermions |2^, |2^, |3(|. The particle 
and hole states in the system of fermions with the harmonic potential interaction corre- 
spond to the giant gravitons of a 1/2 BPS sector of N = 4 SYM |, §, |, f|. Therefore, 
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our solutions correspond to the coherent states of the matrix model with the harmonic 
potential, i.e. to the quasi-classical CFT duals of the giant gravitons in AdS constructed 
in Ref . || . The nonexistence of the quasi-classical CFT dual of the single giant graviton on 
the sphere <S 5 is reflected throught the relation ( |2.20| ) , from which it follows that the soliton 
with M = 1 in the A > 1 case is non-normalizable since M must exceed N, in accordance 
with the conclusion of Ref. § . 

At this point we would like to emphasize a simple relation between systems with A < 1 
and those with A > 1. By substituting Xp(x) = a — m(x) into Eq. (|2.7|) for A > 1 (without 
the term Vcot) and by inserting explicit forms of the solutions for the term p H / p, we find 
that the field m satisfies Eq.(2.7) for A' = 1/A < 1 (with the term Vcot). This agrees with 
the result obtained in Ref. [j3l]l in the k-space (p^ — > — m^/A). The difference is that our 



relation is valid for the BPS equations, whereas in Ref. |31[| the duality relation connects 
quantum Hamiltonians. Recall that the solution of the form similar to that given in the 
third row of Table |] was found as a solution of the dynamic equations of motion of the 
Calogero model in Ref. 1 25]. This signalizes that there is a generalization of our method 
for dynamic equations of motion (some progress has been made in Ref.[16|). Finally, we 
stress that the construction of the Riccati equation is possible for the model trapped in the 
harmonic well, providing a new method for analysing this model. 

2.3 Quantum excitations around semiclassical solutions 

To get an insight into the dynamics of quantum excitations, we expand the Hamiltonian 
fll.6j ) around the semiclassical solution 

p{x,t)=p s (x) + d x r]{x,t), (2.21) 

where rj is a small density quantum fluctuation around the soliton solution p s of Eq. ( |2.7| ) . 
The quadratic part of the Hamiltonian can be written in the following form: 



F (2) = i J dx Ps {x)A\x)A{2 



where we have introduced the operators A 

(A-l 



.4 



-7T„ + I 



o x vrAa z r/ 

2 p s 



satisfying the following equal-time commutation relation: 



A(x),A^y) =(1-A)a 



2 S(x-y) 



+ 2AcV 



xy p s {x) x-y 
Using the equation of motion A(x,t) = i[H,A(x,t)], we obtain the equation 



-id t + 



A-l d x p t 
2 p s 



(fi a A) = -\irp s d x (p s A) H . 



Taking the Hilbert transform of this equation, and using Eq.(p.7|), we find 



-idt + 



A-l d x p s 



r 



A-l 



0: 



{ Ps A) H = \-K Ps d x ( Ps A). 



(2.22) 



(2.23) 



(2.24) 



(2.25) 



(2.26) 
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Defining the fields 

*± = pf ± ip„ ^ = G»^) H ± i(p a A), (2.27) 



we find that 4^ satisfies 
idt 



d x + ^-±d 2 A<P ± = 0. (2.28) 



A^ + ^Pcot(f + v 

This equation can be obtained from the Riccati equation ( |2.10| ) by adding 2i^0 ± / (A — 1) on 
the LHS and expanding <fr(x,t) around the solution & s (x), $> (x,t) = $g (x) + d x cj) (x,t), 
keeping only terms linear in (j). Therefore, one can interpret the field as a fluctuation 
around the conformal field <j? s . Solving Eq.( 2.28|) for the solutions found in section 3, we 
obtain the following results: 
-the operator A is given by 

A = tY1 eiUJnt fnA%) [#KK,s + e(-u n )al s ] , (2.29) 

n,s 

where the operators a ntS satisfy 

[an,a, a m , s ') = \u n \L/iv5 nm 8 ss > (2.30) 
and the functions f UtS are orthonormalized with respect to the measure p s (x) as follows: 

L/2 L 

dxp s (x)f n (x)f miS '(x) = —5 nm S ss >] (2.31) 

-L/2 27r 

-the Hamiltonian up to quadratic terms is given by 

H = E + | J^4,»an,« + J2 e (- UJ n)\un\- (2.32) 

n,s n,s 

The functions f n , s and the eigenvalues u; n for all solutions found in section 3 are given in 
Table g where k = \\Kp /(\ - 1)|. 

Comparing the known dispersion laws (|1.2j), ( p..3j ) for the quantum excitations on the 
uniform background, obtained from the Young diagrams in Ref.|| with the dispersion laws 
given m Table | we find an agreement. At the semiclassical level we obtain that the system 
can be in the fluid phase (uniform density) as well as in the crystal-like configuration. So, 
it would be interesting to calculate the phase transition amplitudes and also to calculate 
the correlation functions determined in this approach by the quadratic Hamiltonian. 
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A<1 



g(l-A)coth(t/2) l-cos(q:r+2y>) 

2tt\ cosh t— cos(qx+2ip) 

1-A a 2 
A-7rfe x 2 +b 2 

PO 



X(k +q)(k n +q) 
4(l-A)fc fc n (2fc + g) 



s ±i(qa: + 2(y5) 



1- 



/n,± 

^i(qx+2<p) \ e ±i(k„-k )x 



kQ + q 



\/ 2fc (l-A) ( 1± fc„ < 



1-005(9:1+2^) fc "> fc 
fe o)^ fc n >fc 



o±i(fc?7 



fc(A-l) S i n ht 
27rA cosh t— cos fcx 
A-l 6 

A-7T ^2_(_^2 
PO 



^"2^ (k„+ko +qj (k„ - ho) 

— — (fe^-fey) 



A>1 



2fc (X-l)(l-e-^) 



/2irp<) 



oiitfen + feo) 1 



fe„>— feo 



2 ( fc 

2 K ™ 
2 (^o — ^n) 



Table 2: Excitations around BPS solutions 



3. Duality-based matrix model 



3.1 Semiclassical solutions 



In attempt to go beyond the case of free fermions, in this section we analyse a generalization 
of the hermitian matrix model, introduced in Refs.[24, ||. Actually, this model was first 
formulated in Ref.|32| as a duality-based generalization of the Calogero model, and is 
defined by the Hamiltonian 



A 

K 



M 



M r> 2 1 



N,M , 

2 ^ (xi 
k 2 /X (k 2 /X - 1) 



a=l 



(Z 



(3.1) 



In Ref.Q, it was shown that for A = 1/2 this model arises from the decomposition of the 
hermitian matrix into the sum of symmetric and antisymmetric matrix. Transformation 
into the hydrodynamic formulation, in terms of the density operators p(x), m(x) and 
the corresponding conjugate operators for k = 1, results in the hermitian collective-field 
Hamiltonian j|] 

H=^j dxp(x)(d x ir p (x)) 2 + ^ J dxm(x)(d x ir m {x)) 2 + 



+ Jdx 

+ dx 



2 

mix) 



2X 



A - 1 d x p(x) 

2 p{x) 
1 — A d x m{x 



-jdxp(x)d : 



2 m(x) 
P 



+ 



+ 



J x-y 

Jdy^l + 
J x-y 



f 



m{y) 
x-y 

My) 
i 

x-y 



y 



- I dxm(x)d x 



x=y 



P 



U x=y 



(3.2) 



The terms in the last line in (3.2) are singular counter terms, which do not contribute in 
the leading order in N and M. This model was analysed in Ref.[33], where some solutions 
were constructed. Here, we use the Riccati equation to construct some new solutions in 
this duality-based model. We are looking for the solutions of two coupled Bogomol'nyi 
equations, 



(A - l)d x p - 2irp(\p H + m H ) = 0, 
(1 - X)d x m - 2-Km{\p H + m H ) = 0. 

Based on the duality, we make an ansatz m H = —Xap H /p. Equation 

(A - l)d x p - 2Xixpp H + 2Xairp H = 0. 



Following the method from section 2, we construct the field <I> = p 
the following Riccati equation: 



H 



d x $ 



Xtt 
A^l 



2A7TO T 

I- <I> 

A- 1 



A7rpo 
A-l 



(po - 2a). 



(3.3) 
(3.4) 

now becomes 

(3.5) 

+ ip which satisfies 
(3.6) 
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The general solution of this equation constructed from the constant solution $ = ipo is 

. A-l iqce^ 2X7r(a-p ) 

<f>(x)=ipo r — — —,q = ; \ >0. 3.7) 

att 1 + ce l i x 1 — A 

The solutions for p and m (c = e lr P- u - v ^ | c | < 1) are 

cosh(u -v) + cos{qx + <j>) c , 

p(x) = a — r -. — , mix) = —-^r, (3.8) 

cosh(u + v) + cos{qx + (p) p[x) 

4A7ra sinh u sinh v c sinh(u — v 



1 — A sinh(u + v) ' Xa 2 sinh(n + v] 



, u > v > 0. 



Taking into account the normalization conditions and the compact support [-L/2, L/2], 
q = 2nn/L, we find the following relations: 

N XN + (1-X)n , XN-M 

Po = -r, ol= — , cothn = 2 + — , 

L XL (1 — X)n 

M c M — (1 — X)n , XN + M 

m = — , - = , cothw = — — . (3.9) 

La L (1 — X)n 



From the solution (3.8) taking <fi = ir, sinh(u/2 — v/2) = aq/2, sinh(u/2 + v/2) = bq/2, 



b > 0, and taking the limit q — > 0, we obtain the one-soliton solution 

x 2 + a 2 Xa 2 a 2 ,2 ^ — I7 

p(x) = a , m(x) = Q =0 + "r 0. (3.10) 

+ b z bp(x) Xira 

Here we would like to emphasize that the properties of the solutions of the dual model 
fl3.2j ) resemble the properties of the dual giant gravitons on AdS (A < 1) and on the sphere 
(A > 1) from Refs.0, |, §. 

3.2 Semiclassical solutions in the singular limit 

In this subsection we discuss the existence of the singular solutions of the duality-based 
model and the methods for finding them. Inspired by the one-soliton solution of the model 



introduced in [24|, we proposed the existence of the multi-soliton solutions in the same 
paper. In the paper |33| we discussed these solutions in the context of the new matrix- 
model interpretation of the duality-based Hamiltonian. We started from an assumption 
that there exists the finite form of a multi-soliton solution. In the case of the one-soliton 
solution, the finite form of the solution was known and we noticed that, in the singular 
limit, the conditions on the parameters of the solution were reduced. This indicates that in 
the singular limit, the Bogomol'nyi equations are less sensitive to the details of the precise 
form of the finite solution. In the hope of getting a hint about the finite form of a multi- 
soliton solution, we used a simple ansatz which in itself was not a solution in its finite 
form, but had the same singular limit as the finite solution. We found that the ansatz 
solved the equations in the singular limit. Although the calculation of the contribution 
to the Hamiltonian leads to an ambiguity in order of taking limits, this problem can be 
avoided by introduction of the finite support. Now, we would like to confirm the existence 
of singular solutions and the correctness of the calculations performed in the paper [33] 
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by using the finite form of the solution constructed in subsection 3.1. Taking the limit 
u — v = 2e — ► of the solution ( |3.8| ) , we find 



( ^ cos 2 (^) (i-A)g 

p(x) = a -7— , a = cotn v, 

sinh 2 ?; + cos 2 (^) 2Avr 

00 i'Zi -\- \)tT (j) 

mix) = (1 — A) y Six — xA, Xi = . 

. q 

i=—oo 

Using the product representation for the trigonometric function we find that the oo-soliton 
solution exactly matches the form proposed in our paper [33]. Furthermore, by taking the 



compact support of lenght L, we find the finite number of solitons. In the case that q 
is small, we find that this solution reduces to the rational ansatz, as proposed in papers 



|2J, |3|. 



In the comment made by V. Bardek and S. Meljanac in [34| regarding this part of our 



published work [33], the authors claim that such solutions do not exist. The authors based 
their claims on the symbolic manipulation with delta functions assuming the independece 
of the regularization procedure. As we are aware that solving nonlinear coupled equations 
is sensitive to the regularization used, the application of the symbolic manipulation is not 
suitable and does not give a unique answer because the product of distributions is not 
uniquely defined. The authors also claim that their paper |35| gives convicing arguments 
that a singular form of the solution does not exist, but they simply have not checked the 
limit c — > in that paper. 

3.3 Open-closed string duality? 

Solutions of the equations obtained in the collective-field formulation of matrix models are 
interpretated as giant gravitons. Although the existence of multi-soliton solutions in their 
finite form is also an interesting question from the point of view of integrable systems, the 
limiting form when one of the fields behaves as a sum of delta functions deserves special 
attention from the point of view of matrix theory. Namely, in this singular configuration, 
it is reasonable to describe the m(x) part of the Hamiltonian by a discrete set of variables. 
Furthermore, we can interpret the continuum part of the Hamiltonian as closed string 
theory and the discrete part of the Hamiltonian as open string theory [Q]. So, by studying 
the relation between these two Hamiltonians we might get a better insight into the open- 
closed string duality. 

The matrix model ( |l.l| ) (the Calogero model) possesses the symmetry of the action 



generated by the operators closing the su(l, 1) algebra |32], g|: 

[T+,T_] = -2T , [T ,T±]=±T ± , 
The generators of the algebra in the discrete case are 




T ({x};A) =-\(Y / x i d i +E ), T„({x};A) = ij> 2 , (3.11 
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and in the continuous, large-N case, 

r+([p]; A) =~J dxp(x)(d x Tr(x)) 2 +l J dx (( A - l )dxP(x) + 2\p(x)jdy^^d x ir(x), 

T (H;A) =~\UJ dxxp{x)d x n(x) + E j , T_ ([/>]; A) = l -j dxx 2 p(x), (3.12) 

where the ground-state energy is Eq(N, A) = (XN(N — 1) + iV)/2. 

Furthermore, it was shown in |32|, ^J, ^] that there existed a strong-weak coupling 
duality in the Calogero model and here we briefly review the main results. In the following 
we use the abbreviation [•] for the arguments of operators, depending on which case is 
under consideration (discrete or continuous) and analogously for the dual system [o]. We 
also need the following definitions: 

lnJ(H,A)={ A Ei^ln(x J -x j ) discrete 

| (A — 1) J dxp(x) In p(x) + Xff dxdyp(x)ln\x — y\p(y) continuous, 



lnV([.],[°]H 



Y^i a ln(xj — Zq) discrete — discrete 

lim^o J J dxdzp(x) ln(x — z — i£)m{z) continuous — continuous(3.14) 
lim e _»o Yla I dxp(x) ln(x — z a — ie) continuous — discrete. 



The strong- weak duality is displayed as follows: 

Zb(M> X)V([m], [o]) = |-T ([o], 1/A) - l -[NM + E (N, A) + E (M, 1/A)]| V([«], [o]), 
r + ([.], A)V([.], [o]) = -AT + ([o], l/X)V([.], [o]). (3.15) 



The operators T + ([»], A) are equivalent to the Hamiltonians (1.4), ( |1.6| 



H C M(n\) = -JH[*},X)T+([*};\)J-*([*],\) (3-16) 

while To([»],A) are equivalent to the same Hamiltonians with an additional harmonic- well 
potential (Calogero-Sutherland Hamiltonians): 

F cs (H;A,a;) = -2wj2 ([•] , A)5([«] ; A, u)T ([•] ; A)^" 1 ([•] ; A, w) J" ^ ( [•] , A) = 

= if C M(M,A)+w 2 r_([.];A), (3.17) 

where 

S-([.];A,o;) = e -^-(W;A) e -^T + ([. ]; A)_ (3 lg) 

Conseqently, we interpret 

H CM ([o]; 1/A) = - J3([o]; l/A)T + ([o] ; 1/A).T S([o]; 1/A) (3.19) 

and 

ifc 5 ([o];l/A,o;/A)=^ji([o] ) l/A)5([o];l/A,a;/A)ro([o];l/A)5- 1 ([o];l/A,a;/A)j4(^ 

= H CM ([o], 1/A) + ([o]; 1/A) , (3.20) 
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as Hamiltonians of the dual system. In the papers [32, 24, 36 1 the duality relations were 
used for construction of the spectrum generating algebra for the two-family model, while 
here we show that they can be used for construction of the eigenstate of one system from 
the known eigenstate of the dual system. Suppose that we know the eigenstate ip{E} with 
energy E, satisfying 

H CS ([*];\u)^{E}=E^{E}- (3-21) 
Then, the eigenstate </>/ ^ of the dual system, with energy E, 

H CS ( [o] ; 1/A, u/X) ^ = , (3.22) 

is given by 

^ } = ^(M,l/A)e-xr-(W;i/A) J ^ {E} jH[.],X)e-» T -W> x W([.]M, (3-23) 
where the integration is performed over the corresponding configuration space and 

XE + E = uj[NM + E (N, A) + E (M, 1/A)]. (3.24) 
The corresponding dual relation also holds for the Hamiltonians Hcm because the eigen- 



states of these are realized as coherent states of the Hcs [^l; 2£, 30], or in another approach 



by the use of unconventional separation of variables in the Schrodinger problem |p7[. Now, 



the quantum mechanics (1.4) of the eigenvalues Xj of M was regarded as an open string/D- 
brane description of the corresponding string theory ||38|| . On the other hand, the collective 
field defined in the large- N limit of the matrix model (3.19) represents the closed string 



excitations ]39j. Inspired by this interpretation we propose an explicit realization of open- 
closed string duality. The relation ( |3.23D tells us how to construct a wave functional of 
closed string excitations described by the collective field p from the wave function of dual 
M open string excitations, and vice versa. 

4. Conclusion 

Motivated by the relevance of soliton solutions in matrix models for giant graviton physics, 
we have addressed the problem of existence of these solutions and have analysed their 
properties. We have introduced a powerful method for obtaining these BPS solutions in 
the collective-field approach. The method is based on the construction of a boundary 
conformal field out of the density of eigenvalues of a matrix, satisfying the Riccati differ- 
ential equation. This method extends to the related Calogero models, in particular to the 
model with the harmonic-well potential. Furthermore, we have established the relation 
between the hydrodynamic Benjamin-Ono equation and the Riccati equation, suggesting 
that our method could be extended to non-BPS equations by inclusion of dynamics. Such 
extension might shed more light onto the underlying boundary conformal theory, an issue 



indicated in Ref.[16|. The solutions we have obtained using this method are connected 



by the duality relation Xp(x) = a — m(x), where m(x) satisfies the BPS equation for the 
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A' = 1/A case. Owing to the su(l, 1) dynamical symmetry these solutions correspond to 
the quasi-classical CFT duals of the giant gravitons on AdS*, (A < 1) and on the sphere 
S 5 (A > 1) H|, [u], |l|, §]. Further, the method has enabled us to solve the equations which 
govern the dynamics of quantum excitations around the uniform and various x— dependent 
backgrounds. These excitations described by the Hamiltonian in quadratic approximation 
represent quantum corrections to the semiclassical solutions. As an application of the re- 
sults for quantum excitations on various x— dependent backgrounds, one could determine 
the correlation functions, the wave functionals of different states and the transition am- 
plitudes among them. On the other hand, one could evaluate the same quantities using 
random matrix theory |^(| (at least for A = 1, 1/2, 2) and then compare the results. 

We have found the finite form of the n-soliton solution in the duality-based matrix 
model, indicating the complete integrability of this model. Owing to the origin of the 
duality-based matrix model and to the properties of its semiclassical solutions, we interpret 
this matrix model as a model of interacting giant gravitons having cubic interaction. The 



BPS solutions for fields p(x) and m(x) (3J3) related by p(x)m(x) = c admit an interesting 
singular limit p{x)m{x) — * 0; at the places where one field is different from zero the other 
field is vanishing. This "black/white" distribution is a characteristic of the two-dimensional 
droplet model of electron gas [41]. Finally, the singular limit of the n-soliton solution 
has motivated us to propose a realization of open-closed string duality. We have found 
the explicit mapping between the eigenstates of the collective-field Hamitonian and the 
eigenstates of the Hamiltonian describing quantum mechanics of nonrelativistic fermions. 
The proposal could be made more precise by expanding the relation ( |3.23|) between the 
wave functional of the closed string excitations and the wave function of dual open string 
excitations around one of the soliton solutions found. We hope to address this issue in a 
future publication. 
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